/* mv.c

   Matrix-Vector multiplication implementation file

   Basile Clout, September 2007
   */

#include "mv.h"

int mvd(const double **mat, const double *vec, double *res, const int m, const int n){

  int i,j, k;

  for (i=0;i<m;i++){
    res[i]=0;
    for(j=0;j<n;j++)
      res[i] += mat[i][j]*vec[j];
  	for(k=0;k<LOOP;k++);		/* Empty computational loop */
  }
  return 1;
}
  


